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Abstract -I establish the fundamental lower bound on the mass that can be measured with a 
nano-mechanical resonator in a given quantum state based on the fundamental quantum Cramer- 
Rao bound, and identify the quantum states which will allow the largest sensitivity for a given 
maximum energy. I show that with existing carbon nanotube resonators it should be possible in 
principle to measure a thousandth of the mass of an electron, and future improvements might 
allow to reach a regime where one can measure the relativistic change of mass due to absorption 
of a single photon, or the creation of a chemical bond. 
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Introduction. — High-quality nano-mechanical res- 
onators can act as extremely sensitive sensors of adsorbed 
material. Impressive progress has been made in this direc- 
tion over the last few years: In 2004, experiments reached 
a level of sensitivity of femto-grams pQ, atto-grams [2], 
and two years later already zepto-grams [3]. Gas chro- 
matography at the single molecular level was achieved a 
year ago [4] , and brought a vast range of chemical and bi- 
ological applications in reach. A mass sensitivity as small 
as half a gold atom has been demonstrated using a nano- 
mechanical resonator based on a carbon nano-tube [5] . At 
the same time, large efforts have been spent to cool down 
a nano-mechanical resonator to its ground state, with the 
ultimate goal of engineering arbitrary quantum states (see 
e.g. [6UT0]). The ground state was reached very recently 
for a piezo-electrical device [TT|. It is therefore natural to 
ask whether the sensitivity of mass measurements could 
be increased further by engineering the quantum state of 
a nano-mechanical resonator, and what would be the truly 
fundamental lower bound on the mass that can be mea- 
sured based only on the laws of quantum mechanics. Early 
on, theoretical investigations tried to find the limitations 
of mass measurements with a nano-mechanical resonator 
[12HT4] . But the bounds which were derived so far as- 
sume that one measures the linear response of the res- 
onator driven at its resonance frequency |12H15| . In the 
experiments, a variety of different read-out and/or cool- 
ing techniques (e.g. optical [TrjH2"2~| . through electrostatic 



effects [2"l l2"3"H2"5] , mechanical [55] , or even field emission in 
the case of a nano-tube [57]) were used. Most of these do 
use linear transduction, but whether this is the optimal 
measurement procedure is an open question. 

The truly fundamental lowest (but achievable) bound 
on the mass sensitivity is a function of the quantum 
state of the resonator, and optimized over all possible 
measurement procedures. It will be calculated below 
using quantum parameter estimation theory, which 
leads to the ultimate limit of sensitivity, the quantum 
Cramer- Rao bound [28]. It becomes relevant once all 
other limitations such as technical noise, adsorption- 
desorption noise, momentum exchange noise, etc. have 
been eliminated |I5j . I will even assume a harmonic 
oscillator without any dissipation (and thus decoherence 
effects), as mixed states can only decrease the ultimate 
sensitivity compared to the pure states from which they 
are mixed [29]. Nevertheless, the bounds I calculate are 
attainable in principle if the idealized conditions are 
met, and therefore set an important benchmark to which 
the performance of existing sensors should be compared 
to. As a guide to further improving the sensitivity 
of mass-sensing using quantum-engineered states of a 
nano-oscillator, I determine the optimal quantum state 
for a given maximum number of excitation quanta in the 
oscillator. 
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Quantum parameter estimation theory. For 

small enough excitation amplitudes, the nano-mechanical 
resonator can be modelled as a harmonic oscillator with 
mass M and effective spring constant D [5] , resonance fre- 
quency oj = y/D/M, and hamiltonian H u = hjj{al } a bJ + h) 
with the usual raising (lowering) operators aj, (o u ). If a 
small mass SM is added to the oscillator, its frequency 
changes to Cj = u(l — e) with e = (1/2)SM/M, and we ob- 
tain the new hamiltonian Hq from H u by replacing oj — > Cj 
everywhere. An arbitrary initial quantum state po is thus 
propagated to p(oJ,t) = U(uj 7 t)poU^ (u,t) (or p(Q,t), re- 
spectively), if no mass (or the mass 5M) is adsorbed at 
t = 0, where U(oj,t) = exp{—iH u t/h). Note that this as- 
sumes that the energy of the oscillator is conserved in the 
adsorption process, i.e. the additional mass is deposited 
with zero differential speed onto the oscillator. The dis- 
tinguishability of the two states p(uj,t) and p(Cj,t) deter- 
mines the smallest SM that can be measured. In general, 
for any density matrix p(x) that depends on some pa- 
rameter x, the smallest Sx that can be resolved from N 
measurements of an observable A (starting always from 
an identically prepared state) is given by [28] 



Sx 
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It has the interpretation of the uncertainty of A in state 
p{x) as judged by N measurements, renormalized by the 
"speed" by which the mean value of A changes as function 
of x. In other words, x has to change by an amount that 
moves the average value of A by at least its uncertainty. 
Optimizing over all possible measurements leads to the 
quantum Cramer-Rao bound [28], 



Sx > Sx„ 



1 



2VN 



3 (p(x),p(x+dx)) 
dx 



(2) 



where dBuies(p{x), p(x + dx)) is the Bures distance be- 
tween p(x) and p(x + dx) (also called Fisher information), 

defined as (feuresCPi,^) = V^Jl - X /F(pi,p 2 ) through 

the fidelity y/F{ Pl ,p 2 ) = tr((^ /2 p 2 p\ /2 )V 2 ). Thus, in 
our case, we obtain the minimal measurable mass SM m i n 
by evaluating the Bures distance between p(u>,t) and 
p(u),i) in the limit e — > 0. It is important to note that ([2]) 
is, in the limit of large N, an achievable lower bound 28 . 
It does take into account effects such as back-action and 
the quantum noise of the system. 

Pure states. — In the case of two pure states, we 
have simply F , \4>){(j>\) 1 ^ 2 = \(if)\4>)\. Starting from 
an initial state | -0(0)) = S^Lo c «l n )w we nave the overlap 
at time t, 

m)\m) = E <c k e-^~- E ^' h 

n,m : k 

xRuu,{m, k)Rs, u (m,n) , (3) 



where i? 1 i ti; (m,n) = &(m|n) w = R U £j{n,m) denotes the 
overlap matrix element between energy eigenstates of the 
two oscillators with frequency Cj and w, and the coeffi- 
cients c* , Cfc are expressed in the energy eigenbasis of the 
unperturbed oscillator. They are [50] 
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if m, n are both even or both odd (otherwise Ruu(m, n) = 
0), and [m, n] denotes the smaller of the two integers m,n. 
The sum over r runs over even (odd) integers for [m, n] 
even (odd), respectively, and y = (u) — Cj) / (cj + Cj) , q = 
2{ujCj) 1 / 2 /{uj + Cj). We need KV'WIV'C*))! to second order 
in e. We find F = 1 + e 2 f{{c m },t) with 



f({Cm},t) 
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and r = wi. Inserting in deuros, we find immedi- 
ately d Bm M(t))(m\, I I) = e|/({c ro },t)| 1 /2, 

and thus 

1 . (6) 

Eq.(j6]) together with ([5]) constitutes the central result of 
this report which we now explore for particular cases. 

Fock state. For \ip(0)) = \n), we have / = -(l/2)(n 2 + 
n + 1) sin 2 r = /^ ock for all n > 0. The largest absolute 
value is achieved for r = ir/2( mod 2tt), and leads to 



5M n 



M 



1 



1 



l!for,»l. (7) 



Thus, one can measure, at least in principle, arbitrarily 
small masses within the same fixed time interval by in- 
creasing the excitation of the harmonic oscillator. In re- 
ality, of course, non-harmonicitics will start to arise at 
some level of excitation and the present analysis will then 
have to be extended to a more complicated hamiltonian 
[31] . The ground state n = of the harmonic oscilla- 
tor allows to measure a mass which, for a single readout, 
can be of the order of the mass of the oscillator itself, 
SM/M > y/2/N/\ sinr|. Increasing r does not help be- 
yond r = 7r/2, as /^ ock is periodic in r. 
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Little Schrddinger cat states. Given that ([5]) depends 
on coherences between states \n), \n + 2), and |n + 4), one 
might wonder whether the precision could be increased 
further by using superpositions of these states. The state 
|V(0)) = (|n) + \n + 2))/y/2 = \ipsx) leads to 

/ = ^-((n + l)(n + 2)sm 2 (2T)-8(n 2 + 3n + 4)sm 2 t)-t 2 . 
16 

(8) 

The maximum of the periodic term is again achieved for 
r = it/2 ( mod 2vr). For fixed r ^ kw (k £ N) and n r, 
this term dominates and leads to <5M m j n /M ~ y/2/N/n, 
just as for the Fock state. However, for fixed n, we can 
get an arbitrarily small SM/M by increasing r, as the last 
term in © leads to 5M min /M ~ xJ2/N/t for r > n. 
Note that this improvement is beyond the usual factor 
1/Vt from increasing the measurement time. Indeed, the 
sensitivity SMy/i in (g/v / Hz) still improves as cx 1/y/t. 
The state |V(0)) = (\n) + \n + 4»/>/2 = |V>S2> gives 

/ = -^(2(" 2 + 5n + 11) sin 2 r (9) 
+ y/(n + l)(n + 2){n + 3)(ra + 4) cos(2r) sin 2 t) 
-4t 2 . 

The maximum of the periodic terms (relevant for fixed r 
and n > r > 1) is close to r = tt/3 with |/| ~ 9n 2 /32, 
which gives a SM m [ n 33% larger than for a Fock state with 
n excitations. For fixed n and r»n, the factor 4 in front 
of t 2 in (|9]) reduces SM m - m /M by a factor 2 compared to 
IVsi)- 

Coherent state. For a coherent state \a) = 
exp(— |a| 2 /2) J2^=o( an / V^-)\ n ) ' an d a £R, we have 

/ = - Q + a 2 ^j sin 2 r - a 2 r (r + sin(2r)) . (10) 

For fixed r ^ 1 and a 2 1, we find / ~ — a 2 r 2 , and 
hence SM mia /M ~ l/(VNorr) = l/( v /iV(n)r). Thus, 
for a coherent state, the sensitivity scales as the inverse 
square root of the (average) number of excitations in 
the oscillator, and one can, at least in principle, resolve 
arbitrarily small masses. We see again that <5M m i n reduces 
faster than \j\ft with measurement time. Compared to a 
Fock state with n = a 2 , there is a factor \J (n)/2 penalty 
in the scaling with (n) (but one gains a factor 1/r). For 
a = we find / = — sin 2 r/2 which leads back to the 
result for the Fock state n = 0. " 

Optimal state. What is the best sensitivity that can 
be achieved for a given maximum number L of excita- 
tions and fixed measurement time? From (JS]) we see that 
for t 3> 1, the terms quadratic in t dominate and give 
simply |/| = T 2 ((n 2 ) — (n) 2 ). Hence, in this case the 
optimal pure state is the one which maximizes the ex- 
citation number fluctuations. One easily shows that this 
state has the form of an "ON" state (half a "NOON" state 
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Fig. 1: Inverse minimal measurable mass M/<5M m i n (for iV = 1) 
as function of r for selected pure states. Green line: Fock state 
| ra = 3); Full black line: ON state |^on) = (|0) + |3))/V2; 
Dashed black line: asymptotic behavior 3r/2 for I^on); Red 
line: optimal state with at most L = 3 quanta in the oscilla- 
tor; Dashed blue line: Coherent state with the same average 
number of excitations as the ON state ((n) = 3/2). 

ED), IV>on> = (|0) + c^\L))/V2, where ip is an irrele- 
vant phase which we will choose equal zero. It leads to 
| /on | = \ f(TpON,T~)\ ~ t 2 L 2 /4, and thus a minimal mass 
(5M m j n /M = 2/( s/NtL). FigQ] shows a comparison of the 
(inverse) minimal mass for "0ON with the true minimal 
mass for given r and the same L, obtained by numerically 
maximizing |/|, for L = 3. We see that |/on| approxi- 
mates the best possible |/| very well, even for r ~ 1. For 
t = kir, k £ N, | /on | gives in fact the exact result, as is 
obvious from ([5]). FigQ] also shows the result for a coherent 
state with the same average number of excitations as the 
ON state, (n) — 3/2. It leads for small r to comparable 
sensitivity as the optimal state with L = 3. Atr = 7r/2, 
the optimal pure state with L = 3 allows still a reduction 
of (5M m i n /Af by ~ 4% compared to Von, and by ~ 18% 
compared to the Fock state with the same L. 

Fig(2]shows the Wigncr function, defined for a pure state 
IV) by [33] 

i r°° 

W(x,p) = - / dyV* (x-y)ip(x + y)e" 2iOT , (10) 

(with all lengths in units of the oscillator length 
x = h 1 / 2 /(DM) 1 /' 1 and p in units of h/x Q ), 
for the optimal state for L = 4 and t = tt/2, 
IVopt) ^ (-0.62057+0.0305i)|0)-(0.00059+0.01198i)|2) + 
(0.78252 - 0.038852i)|4). |Von) and |V' op t) have very 
similar Wigner functions, characterized by four lobes 
in azimuthal direction which guarantee minimal phase 
uncertainty, as is to be expected from the requirement of 
minimal noise and maximum uncertainty in the number 
of excitations. Rotations of |Vopt) through evolution 
with the unperturbed hamiltonian before the adsorption 
of mass clearly leave <5Af m ; n invariant. It is an open 
question how to experimentally realize the optimal state. 
In [11] , non-trivial quantum states were obtained by 
coupling the harmonic oscillator to a super-conducting 
qubit whose quantum state can be readily manipulated. 
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Fig. 2: Wigner function of the optimal initial state for L = 4, 
T = tt/2 

Other schemes for engineering the quantum state of a 
nano-mcchanical harmonic oscillator were proposed in 

[S3- 

Among the pure states considered, the coherent states 
certainly come closest to the typical experimental situa- 
tion, where the oscillator is cooled to low temperature and 
driven on resonance. Inserting typical numbers for micro- 
machined resonators, M = 10~ 16 g [14] . u =lGHz, an evo- 
lution time r = 10 6 , and an excitation with (n) ~ 10 10 
quanta (driving energy Ed = 10~ 15 J in [Hj), we find 
SM m - ln ~ 10~ 27 g/y/N, or roughly the mass of an elec- 
tron for a single readout, N = 1. Higher masses (e.g. 
M ~ 10~ 14 <7 in [3]) give proportionally higher 6M m i n , ev- 
erything else equal. In assuming r = 10 6 , we have made a 
pessimistic estimate in the sense of using the shortest sens- 
ing time allowed by the inverse bandwidth 1 kHz. Measur- 
ing during longer times decreases SM m { n further. On the 
other hand, t might be limited by decoherence and dissi- 
pation, such that for later times another (mixed) quantum 
state becomes relevant. These effects are, at least in prin- 
ciple, avoidable, and in such a highly idealized situation 
(which is, however, relevant as utlimate achievable goal), 
r in eq.([5]) is given by the measurement time. The mass 
SM m i n ~ 10 _27 g agrees with the prediction of ffj], but 
the agreement appears to be a coincidence: The result in 
[14], based purely on noise considerations, still decreases 
as l/\fQ with the quality Q of the resonator, whereas (fT0| 
is independent of Q, taken as infinity in the present analy- 
sis. Also, while in the regime r»a relevant for the above 
numbers (a = 10 5 ) both 8M m - m /M and the result in [14] 
scale as 1/ \J (ra), QT] predicts a proportionality to 1/ y/r if 
one identifies the inverse bandwidth 1/A/ with t, instead 
of the 1/r behavior that follows from ([TO]) . 

Carbon nanotube resonators have typically much 
smaller masses than micro-engineered ones (of order M ~ 
10~ 18 g [5]) with comparable resonance frequency (w = 
2tt x 328.5MHz in [5]), and can therefore resolve in prin- 
ciple even smaller masses. Assuming a coherent state 
with oscillation amplitude of about lOnm for the carbon 
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Fig. 3: Inverse minimal measurable mass M/SMmm in a 
thermal state as function of dimensionless time r = tot 
for inverse dimensionless temperatures z = fa/(feT) = 
0.2,0.5,1.0,2.0,5.0 and 10.0 (red, orange, yellow, green, blue, 
and dashed purple line, respectively). Inset: M/<5M m j n as func- 
tion of 2 for t — tt/2. 

nanotube resonator in [5] and a sampling time of 100ms, 
8M nlin according to (JTUJ) is of the order of a thousandth of 
an electron mass. 

Mixed states. — In general, due to the joint- 
convexity of the Bures distance, mixed states do not allow 
better sensitivities than the pure states from which they 
are mixed, as long as the weights are independent of the 
parameter to be measured and the evolution is linear in 
the state [29] . Evaluating the Bures distance between two 
mixed states is much more difficult than for pure states. 
Nevertheless, one can evaluate the Bures distance numer- 
ically. One may also obtain an upper bound on dsures 
using the joint-convexity of dsuros, which leads to a (typ- 
ically non-achievable) lower bound on <5Af m i n . An achiev- 
able upper bound on 5M m - m can be found by considering 
a particular measurement Am JT]) . 

As an example, consider the thermal state p = 
Lr=oP«l n )H> with P" = e~ nz (l - e" z ), z = (3hoj, 
P = 1/ksT the inverse temperature, and Ub the Boltz- 
mann constant. Using the invariance of the thermal 
state under the time evolution governed by H u and the 
joint convexity of e?Bures(/9i> P2), one shows easily that 

dBurcs/dx < EnPn\fn° Ck \ 1/2 < I Sin r|/( ^2(1 - CXp(-z)). 

For z — > oo, this bound coincides with the exact re- 
sult for the groundstate n = 0. We may choose a 
measurement of the width A = x 2 as a way of mea- 
suring the change of mass. Thermal average (a; 2 ) = 
(H/(2VDM)) cot h(g/2) and fluctuations (S(x 2 )) = {(x 4 )- 
(x 2 ) 2 ) 1/2 = (H/V2DM) cot h(z/2 ) give an achievable up- 
per bound SM min /M < lyjl JN sinh z / (sinh z - z). For 
z — > oo this bound is only a factor 2 above the best pos- 
sible value for the groundstate n = 0, whereas for z — > 
the bound diverges. Fig[3] shows the exact M/SM m i n ob- 
tained by evaluating the Bures distance numerically. We 
see that 5M m - m /M is periodic in r, just as for the ground 
state. Increasing the temperature helps, as higher Fock 
states start to contribute, but at most a factor v2 can be 
gained, and SM mm remains bounded by the mass of the 
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oscillator itself for all temperatures. 

Conclusions. — In summary, I have calculated the 
smallest measurable adsorbed mass 5M m i n on a nano- 
mechanical harmonic resonator in an arbitrary pure state, 
based on the fundamental quantum Cramer-Rao bound. 
The analysis shows that a coherent state allows to achieve 
a SM m i n /M that scales for fixed measurement time as 
the inverse square root of the average number of excita- 
tions. For a given maximum number n of excitations, 
I found the optimal quantum state, which for large t is 
IV'On) = (|0) + \n))/^/2, with a sensitivity that scales as 
1/n. For r ^ kn, k e N, the sensitivity can be further en- 
hanced. Even with a coherent state of a carbon nanotubc 
resonator [5] , the smallest resolvable mass should be of the 
order of a thousandth of an electron mass. If two more or- 
ders of magnitude could be gained (say by increasing r and 
the number of excitations), the regime could be reached 
where one can weigh the relativistic mass change due to 
the formation of a chemical bond or the absorption of a 
photon (energies of order 1 cV) . 
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